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SYSTEMS AND METHODS FOR FILTERING 
IMAGES 

BACKGROUND OF THE INVENTION 

[0001] This invention relates generally to imaging systems and more 
particularly, to systems and methods for filtering images. 

[0002] For computed tomography (CT) applications, a high contrast- 
to-noise ratio is desirable to detect low contrast lesions and high density organs, such 
as bones and stents in hearts of patients. However, boosting the contrast increases the 
noise of a CT image. 

BRIEF DESCRIPTION OF THE INVENTION 

[0003] In one aspect, a method for filtering images is provided. The 
method includes obtaining an image, and obtaining a final pixel value by performing a 
filtering operation on an initial pixel value of at least one pixel of the image and by 
modulating the filtering operation with a gain factor that is a fimction of the initial 
pixel value. 

[0004] In another aspect, a method for filtering images is provided. 
The method includes obtaining a computed tomography (CT) image, and obtaining a 
final pixel value by performing a filtering operation on an initial pixel value of at least 
one pixel of the CT image and by modulating the filtering operation with a gain factor 
that is a fimction of the initial pixel value. 

[0005] In yet another aspect, a computer-readable medium encoded 
with a program is provided. The program is configured to obtain an image, and obtain 
a final pixel value by performing a filtering operation on an initial pixel value of at 
least one pixel of the image and by modulating the filtering operation with a gain 
factor that is a fimction of the initial pixel value. 

[0006] In yet another aspect, a computer is provided. The computer is 

configured to obtain an image, and obtain a final pixel value by performing a filtering 

operation on an initial pixel value of at least one pixel of the image and by modulating 

the filtering operation with a gain factor that is a fimction of the initial pixel value. 
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[0007] In another aspect, a computed tomographic (CT) imaging 
system for filtering CT images is provided. The imaging system includes a detector 
array having a plurality of detectors, an x-ray source positioned to emit x-rays toward 
the detector array, and a processor operationally coupled to the detector array. The 
processor is configured to obtain an image, and obtain a final pixel value by 
performing a filtering operation on an initial pixel value of at least one pixel of the 
image and by modulating the filtering operation with a gain factor that is a function of 
the initial pixel value. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0008] Figure 1 is a embodiment of a CT imaging system in which 
systems and methods for filtering images is implemented. 

[0009] Figure 2 is a block schematic diagram of the CT imaging 
system of Figure 1. 

[0010] Figure 3 is a flowchart of an embodiment of a method for 
filtering images. 

[001 1] Figure 4 is a plot of gain factor versus relative pixel values to 
obtain the gain factor from the relative pixel values. 

[0012] Figure 5 shows CT images that are generated with and 
without using the method of Figure 3. 

[0013] Figure 6 shows CT images that are generated with and 
without using the method of Figure 3. 

DETAILED DESCRIPTION OF THE INVENTION 

[0014] In some known CT imaging system configurations, an x-ray 

source projects a fan-shaped beam which is coUimated to lie within an X-Y plane of a 

Cartesian coordinate system and generally referred to as an "imaging plane". The x- 

ray beam passes through an object being imaged, such as a patient. The beam, after 

being attenuated by the object, impinges upon an array of radiation detectors. The 

intensity of the attenuated radiation beam received at the detector array is dependent 

upon the attenuation of an x-ray beam by the object. Each detector element of the 
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array produces a separate electrical signal that is a measurement of the beam 
attenuation at the detector location. The attenuation measurements from all the 
detectors are acquired separately to produce a transmission profile. 

[0015] In third generation CT systems, the x-ray source and the 
detector array are rotated with a gantry within the imaging plane and around the object 
to be imaged such that the angle at which the x-ray beam intersects the object 
constantly changes. A group of x-ray attenuation measurements, i.e., projection data, 
from the detector array at one gantry angle is referred to as a 'View". A "scan" of the 
object comprises a set of views made at different gantry angles, or view angles, during 
one revolution of the x-ray source and detector. 

[0016] In an axial scan, the projection data is processed to constract 
an image that corresponds to a two dimensional slice taken through the object. One 
method for reconstructing an image from a set of projection data is referred to in the 
art as the filtered back projection technique. This process converts the attenuation 
measurements from a scan into integers called "CT numbers" or "Hounsfield units", 
which are used to control the brightness of a corresponding pixel on a cathode ray 
tube display. 

[0017] To reduce the total scan time, a "helical" scan may be 
performed. To perform a "helical" scan, the object is moved while the data for the 
prescribed number of slices is acquired. Such a system generates a single helix from a 
one fan beam helical scan. The helix mapped out by the fan beam yields projection 
data from which images in each prescribed slice may be reconstructed. 

[0018] Reconstruction algorithms for helical scanning typically use 
helical weighing algorithms that weight the collected data as a fimction of view angle 
and detector channel index. Specifically, prior to a filtered backprojection process, the 
data is weighted according to a helical weighing factor, which is a fiinction of both the 
gantry angle and detector angle. The helical weighting algorithms also scale the data 
according to a scaling factor, which is a function of the distance between the x-ray 
source and the object. The weighted and scaled data is then processed to generate CT 
numbers and to construct an image that corresponds to a two dimensional slice taken 
through the object. 



[0019] As used herein, an element or step recited in the singular and 
preceded with the word "a" or "an" should be understood as not excluding plural said 
elements or steps, unless such exclusion is explicitly recited. Furthermore, references 
to "one embodiment" of the present invention are not intended to be interpreted as 
excluding the existence of additional embodiments that also incorporate the recited 
features. 

[0020] Also as used herein, the phrase "reconstructing an image" is 
not intended to exclude embodiments of the present invention in which data 
representing an image is generated but a viewable image is not. However, many 
embodiments generate (or are configured to generate) at least one viewable image. 

[0021] Referring to Figures 1 and 2, a multi-slice scanning imaging 
system, for example, a computed tomography (CT) imaging system 10, is shown as 
including a gantry 12 representative of a "third generation" CT imaging system. 
Gantry 12 has an x-ray source 14 that projects a beam of x-rays 16 toward a detector 
array 18 on the opposite side of gantry 12. Detector array 18 is formed by a plurality 
of detector rows (not shown) including a plurality of detector elements 20 which 
together sense the projected x-rays that pass through an object, such as a medical 
patient 22. Each detector element 20 produces an electrical signal that represents the 
intensity of an impinging x-ray beam and hence the attenuation of the beam as it 
passes through object or patient 22. During a scan to acquire x-ray projection data, 
gantry 12 and the components moimted thereon rotate about a center of rotation 24. 
Figure 2 shows only a single row of detector elements 20 (i.e., a detector row). 
However, multislice detector array 1 8 includes a plurality of parallel detector rows of 
detector elements 20 such that projection data corresponding to a plurality of quasi- 
parallel or parallel slices can be acquired simultaneously during a scan. 

[0022] Rotation of gantry 12 and the operation of x-ray source 14 are 

govemed by a control mechanism 26 of CT system 10. Control mechanism 26 

includes an x-ray controller 28 that provides power and timing signals to x-ray source 

14 and a gantry motor controller 30 that controls the rotational speed and position of 

gantry 12. A data acquisition system (DAS) 32 in control mechanism 26 samples 

analog data from detector elements 20 and converts the data to digital signals for 

subsequent processing. An image reconstructor 34 receives sampled and digitized x- 

ray data from DAS 32 and performs high-speed image reconstruction. The 
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reconstructed image is applied as an input to a computer 36 which stores the image in 
a mass storage device 38. 

[0023] Computer 36 also receives commands and scanning 
parameters from an operator via console 40 that has a keyboard. An associated 
cathode ray tube display 42 allows the operator to observe the reconstructed image 
and other data from computer 36. The operator supplied commands and parameters 
are used by computer 36 to provide control signals and information to DAS 32, x-ray 
controller 28 and gantry motor controller 30. Li addition, computer 36 operates a table 
motor controller 44 which controls a motorized table 46 to position patient 22 in 
gantry 12. Particularly, table 46 moves portions of patient 22 through gantry opening 
48. 

[0024] In one embodiment, computer 36 includes a device 50, for 
example, a floppy disk drive, CD-ROM drive, DVD drive, magnetic optical disk 
(MOD) device, or any other digital device including a network connecting device such 
as an Ethernet device for reading instructions and/or data from a computer-readable 
medium 52, such as a floppy disk, a CD-ROM, a DVD or an other digital source such 
as a network or the Internet, as well as yet to be developed digital means. In another 
embodiment, computer 36 executes instmctions stored in firmware (not shown). 
Computer 36 is programmed to perform functions described herein, and as used 
herein, the term computer is not limited to just those integrated circuits referred to in 
the art as computers, but broadly refers to computers, processors, microcontrollers, 
microcomputers, programmable logic controllers, application specific integrated 
circuits, and other programmable circuits, and these terms are used interchangeably 
herein. 

[0025] Figure 3 is a flowchart of an embodiment of a method for 

filtering a CT image. The method is executed by computer 36 after it receives data 

representative of a CT image from image reconstructor 34. The method includes 

obtaining 60 a CT image and determining 62 a threshold value J. The method fiirther 

includes generating 64 a gain factor curve, which is shown in Figure 4, as a fimction 

of a relative pixel value P/i, j) of a pixel (i, J) of the CT image. A gain factor has 

positive and negative values, such as, for instance, ranging from -0.5 to 1. The gain 

factor curve varies with CT applications. For example, stent enhancement in a cardiac 

application generates a different gain factor curve than a gain factor curve generated 
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in an inner acoustic channel application. The method further includes calculating 66 
an effective pixel value P/ij) from a pixel value P(iyj) of the pixel (ijj) by using 

Pe(iJMP(iJ)^P(i-lJ)+P(i^hj)+P(iJ'l)^^^^ 

Equation(l) 

where P(i-h jh Pft-^h P(h j-V, and Pfi, j+1) are pixel values of 
pixels that are within 1 xmit of the pixel (7, j). In an alternative embodiment, Pe(W is 
obtained from pixel values of pixels that are within n units of the pixel (Uj), where n 
is a positive integer. Step 66 reduces the impact of noise for the pixel (7, j) of the CT 
image by the averaging operation. Step 66 is performed for every pixel of the CT 

image so that noise in each of the pixels of the CT image is reduced. 

[0026] The method further includes calculating 68 a relative pixel 
value P/ij) from the effective pixel value by using 

Pr(hj) =Pe(hj)/T Equation(2) 

P/ij) is limited to being equal to or less than 1. The method further 
includes calculating 70 the gain factor for the pixel (iyj) by using the gain factor curve 
of Figure 4. In an embodiment, the gain factor is obtained as 

Gain (i, /) = "0.35 + 0.1 * (i, y )+ 0.15 * (i, jf + 0.2 * (i, jf + 0.4 * (i, jY + 0.5 * P, (i, jY 

Equation(3) 

where * is a multiplication operation. 

[0027] The method further includes categorizing 72 the CT image 
into at least two regions of low, medium, and high density. The CT image is 
categorized based on a pixel where a gain factor, as shown in the gain factor curve, is 
zero. The number of regions depends on the CT application. The method further 
includes obtaining 74 a final pixel value Pf{ifj) of the pixel (/,y) by using 

Pf ihj) = P{ijh{Piijh decon{p{ij))YGain{ij) 

Equation(4) 
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where decon{p{fj) ) is a deconvolution or a filtering operation 
performed on the initial pixel value The gain factor modulates the filtering 

operation. For pixels of the CT image whose values are less than a product of a 
constant and the threshold value J, the deconvolution operation is a smoothing 
operation, and for pixels of the CT image whose values are greater than the product, 
the deconvolution operation is a sharpening operation. The constant is determined 
from the gain factor curve. The smoothing operation is limited to pixels that belong 
to a particular region of the at least two regions to maintain a low contrast resolution 
and to avoid oversmoothing over different structures in the CT image. The region in 
which the smoothing operation is applied depends on the CT applications. The 
threshold value T may be adjusted later on the fly by the user based on the CT 
applications to change smoothness and sharpness of the CT image. 

[0028] Figure 5 shows CT images 80 and 82 of stents of patient 22. 
Image 80 is generated without using the method of Figure 3 and image 82 is generated 
by using the method of Figure 3. Figure 6 shows CT images 84 and 86 of an inner 
acoustic canal of patient 22. Image 84 is generated without using the method of Figure 
3 and image 86 is generated by using the method of Figure 3. Images 82 and 86 have 
better contrast than images 80 and 84. 

[0029] Hence, the herein described systems and methods provide a 
generalized non-linear post-processing image filter for CT applications. The filter 
enhances the contrast for high-density objects, such as, bony structures and stents, in 
CT images while reducing image noise for soft tissues. This enables the user of CT 
system 10 to maintain or lower scan techniques and simultaneously obtain higher 
spatial resolution for the high-density objects. By reducing the noise in the soft 
tissues, the user increases the contrast-to-noise ratio, which improves detectability of 
low contrast regions. 

[0030] While the invention has been described in terms of various 
specific embodiments, those skilled in the art will recognize that the invention can be 
practiced with modification within the spirit and scope of the claims. 
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